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The DAM A/LIBRA collaboration has detected an annual modulation of the recoil rate in NaT 
crystals with the phase expected for WIMP scattering events. This signal is dramatically inconsistent 
with upper limits from other experiments for elastically scattering weak-scale WIMPs. However, the 
results are compatible for the case of inelastic dark matter (iDM). The iDM theory, as implemented 
by Tucker-Smith and Weiner, constrains the WIMP to a tight contour in o"„ — (5 space, where 5 is the 
mass difference between the ground state and excited WIMPs. An urgent priority in direct detection 
is to test this scenario. The crucial test of the iDM explanation of DAMA - an experimentum crucis 
- is an experiment with directional sensitivity, which can measure the daily modulation in direction. 
Because the contrast can be 100%, it is a sharper test than the much smaller annual modulation in 
the rate. We estimate the significance of such an experiment as a function of the WIMP mass, cross 
section, background rate, and other parameters. The proposed experiment severely constrains the 
DAMA/iDM scenario even with modest exposure (~ 1000 kg • day) on gaseous xenon. 

PACS numbers: 95.35.-|-d 



I. INTRODUCTION 

The DAMA claim [l| of an annual modulation signal 
has long appeared to be in conflict with non-detections 
in other experiments 0- Though recent limits by 
XENONIO $,^ and CDMS II [3] appear to rule out the 
DAMA region of parameter space by a factor of 100 in 
cross section, DAMA/LIBRA Q has recently confirmed 
their previous annual modulation result and increased 
the significance to 8.2a. This confiict has motivated se- 
rious discussion of models beyond the simplest elastic 
scattering of weak-scale WIMPs, with the hope of ac- 
commodating DAMA as well as the other limits. 

At least four approaches have been considered: 1. elec- 
tron scattering [3]; 2. spin dependent scattering 1C|; 
3. light dark matter [111]; and 4. inelastic scattering 12 1. 
The first hypothesizes that the signal in DAMA is scat- 
tering of WIMPs off of electrons. Significant momentum 
can be transferred to the electron during the small frac- 
tion of the time (< 0.1%) that it finds itself near the 
nucleus and at moderately relativistic speeds. However, 
this small fraction must be balanced by an uncomfort- 
ably large cross section, which is almost certainly ruled 
out by early Universe (CMB) constraints. 

The spin-dependent scattering argument attempts to 
circumvent limits from CDMS in Si for example by posit- 
ing that the cross section is strongly dependent on nuclear 
spin. However, recent experiments [l3| have significantly 
tightened constraints on this scenario, and the allowed 
regions require a significant drop in the background in 
the signal region '14]. While small regions of parameter 
space are still allowed, we do not consider this here. 

Another suggestion is that the DAMA recoil events 
are not in the energy range first suspected. Assuming 
recoils off of iodine, the quenching factor of 0.09 implies 



that the 2 — 6 keVee observed energy corresponds to a 
recoil energy of 22 — 66 keVr. It has recently been sug- 
gested that "channeling" , i.e. alignment of the recoil with 
principal directions in the crystal lattice, creates an ef- 
fective quenching factor of unity for some fraction of the 
events |15|. In this case, there is a small amount of pa- 
rameter space available for lighter WIMPs 5 GcV) 
still compatible with other limits [ij, [IB]- In general, 
light WIMPs have difficulty with constraints from the en- 
ergy spectrum of the unmodulated DAMA signal [13, El]- 
While further exploration of light WIMPs may be war- 
ranted, we do not consider this option here. 



A. The DAMA/iDM Scenario 

The inelastic scattering scenario of Tucker-Smith & 
Weiner [l^, [H, [13] takes a different approach: inelastic 
dark matter (iDM) has an excited state some 6 ~ 100 keV 
above the ground state. The origin of this excited state is 
unimportant for the present arguments; see [2l| for one 
realization of this idea. Elastic scatterings off of nuclei 
are suppressed by at least two orders of magnitude with 
respect to the inelastic scatterings, leading to a preferred 
energy threshold with few events at low energies. The 
high sensitivity of e.g. XENONIO to low-energy scatter- 
ings (which dominate in the standard elastic scattering 
models) means that even a small exposure time (316 kg 
day) can place record-beating limits on the elastic cross 
section. Because iDM does not produce such low-energy 
events, it is plausible that the much larger combined ex- 
posure time of DAMA/LIBRA and DAMA/Nal (300,000 
kg day) could see the higher energy events invisible in the 
other experiments. 

Models of iDM are simple to construct, for instance 



2 



a four th-g eneration neutrino, coupling through the Z- 
boson |19| . a mixed sneutrino KK states in RS the- 
ories [22], in composite models [23], or in theories with 
light mediators 24]. In fact, off-diagonal couplings are 
very natural in dark matter theories, with only the small 
splitting S remaining to be explained. 

In an annual modulation experiment, iDM enjoys an 
additional enhancement relative to elastic models be- 
cause only WIMPs on the high velocity tail scatter. The 
modulation can be much larger than the 2-3% expected 
for elastic scattering, partially compensating for the fact 
that the majority of WIMPs are below threshold and do 
not scatter. 

If the direct detection data from DAMA and others are 
taken at face value as nuclear WIMP scattering events, 
they argue strongly for further experiments designed to 
test iDM. The experiment must make predictions beyond 
the already observed annual modulation so that a posi- 
tive result would add substantially to the believability of 
the result. Such a make-or-break experiment is known 
as a "critical experiment," or experimentum crucis [5l| . 
In the next section we describe such an experiment and 
discuss the limits obtained. 



B. Advantages of Directional Sensitivity 

The DAMA result is compelling enough to motivate 
further experiments involving iodine or other nuclei of 
similar mass. Direct detection experiments generally fall 
into 3 categories, based on their background rejection 
strategy. Some (CDMS II, XENONIO, etc.) reject indi- 
vidual electron scattering events and look for the residual 
signal from WIMP scattering. 

Another strategy for dealing with background is to 
search for the annual modulation of the signal (DAMA) 
brought about by the Earth's velocity around the Sun, 
added to the velocity of the Sun around the Galaxy. 
The assumption is that the WIMP velocities are nearly 
isotropic, and the Sun moves through the WIMPs at 
roughly 200 km/s. The Earth moves around the Sun 
at Vorb ~ 30 km/s in an orbit inclined hy i ~ 60° with 
respect to the Sun's velocity, introducing a modulation 
of Uor& cos(i) ~ 15 km/s. This method has the virtue of 
ignoring all steady state instrumental backgrounds, but 
is vulnerable to backgrounds that vary with the seasons. 
Though DAMA has placed stringent limits on variations 
in temperature, humidity, radon gas, line voltage, and 
anything else known to vary by season Q, this remains 
a persistent concern. 

A third strategy is to use directional information [2^ . 
Because the scattering events should originate, on aver- 
age, from a specific direction on the sky {£ — 90°, 5 = 0°), 
a daily modulation in direction due to the rotation of the 
Earth is a sharp test of the WIMP scattering model. As 
with the annual modulation, many other backgrounds 
may be expected to vary on a daily timescale. However, 
as the Earth orbits around the Sun, the angle between the 



Sun direction and the WIMP signal varies from 60° (^^7 
March) to 120° (~9 September). Also, any Sun- related 
oscillation (365.25 yr~^) is orthogonal to the WIMP sig- 
nal (366.25 yr^^) over one year. This separation allows 
a much sharper test than the annual modulation alone, 
even in the limit of low statistics. Furthermore, direc- 
tional detectors have excellent background rejection and 
can distinguish between recoils of nuclei and other parti- 
cles by correlating the length and energy of recoil tracks. 

In the context of iDM, a directional experiment has 
another advantage. The minimum velocity Vmin for a 
WIMP to scatter with a nuclear recoil of energy Efj is: 

where fj, is the nucleus- WIMP reduced mass 
m-^rriN / (rn-^ + ttin). Because of the energy thresh- 
old, most events result from WIMPs in the tail of the 
WIMP velocity distribution, and therefore most events 
happen near threshold. This is advantageous because 
events at threshold have a sharply peaked angular 
distribution, making the directional discrimination even 
more pronounced. The energy-dependent maximum 
recoil angle is 

COS -f^ax{E li) = (2) 

VE 

Here 7 is the angle between the velocity of the Earth 
and the recoil velocity in the Earth frame, and v^sc is 
the Galactic escape velocity from the Solar neighbor- 
hood. For the benchmark models considered here, 7 is 
constrained to be within '^lOO degrees of the Earth's di- 
rection. Furthermore, as with annual modulation, the 
total number of events should vary through the year in a 
predictable way. These advantages allow a decisive test 
of the DAMA/iDM scenario with modest experimental 
effort. 

In this article, we evaluate the sensitivities for the 
DAMA/iDM scenario as a function of WIMP mass, delta, 
and other parameters. We focus on a set of benchmark 
models, given in TablelH that can simultaneously explain 
DAMA and satisfy constraints from other experiments 
[20| . Note that the = 70 GeV benchmark cannot ac- 
tually explain the DAMA data because of the predicted 
asymmetry in the modulation amplitude during summer 
and winter. However, we include the benchmark as a 
worst-case scenario. These models should give all the 
general features and sensitivities (within an order of mag- 
nitude) for the entire available parameter space of iDM, 
though there is some flexibility in the WIMP parameters 
due to the uncertainty in the halo distribution and as- 
trophysical parameters. We find that in most parts of 
parameter space, 1000 kg days of exposure is sufficient 
to confirm or refute DAMA/iDM at high confidence. 
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FIG. 1: Differential rates dR /{dE ad cos j) for the benchmark models given in TableUfor Vesc ~ 500 km/s, as well as an elastic 
WIMP. In each case, the differential rate is normalized so that the total rate is unity. Outside the region indicated by the 
dashed line, scattering events are kinematically forbidden. 
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TABLE I: Benchmark models for Vesc = 500 km/s, vq = 220 
km/s 20] . In the last row we have listed the benchmark model 

-- 600 km/s. 



II. EXPERIMENTAL SETUP 

Before discussing the specifics of the experiment, we 
can address a few basic questions of exposure and energy 
range. DAMA/LIBRA reports a cumulative modulation 
in the 2 — 6 keVee range of 0.052cpd/kg. The quoted 
energy range is related to the nuclear recoil energy by a 
quenching factor q = Eee/Epfn ~ 0.09 for iodine. Thus, 
2-6 keVee « 22 - 66 keVr. 



In the extreme case where the modulation is 100% (i.e., 
no scattering at all occurs in the winter), the signal is 
essentially directional. One would need approximately 
200 kg • day in the summer to yield 10 events of signal, 
roughly the number of events needed for an unambigu- 
ous detection at zero background, as we will discuss in 
Section ITVAI Consistency with other experiments is also 
possible with ~ 20% modulation [205, with only 40 kg • 
day needed for a clear discovery. 

However, this estimate assumes that the signal occurs 
in an energy range which is detectable at a directional 
experiment, and this, we shall see, is very unlikely to 
be the case. A directional experiment will likely have a 
higher energy threshold. 

The DAMA/LIBRA signal peaks near Er « 3 keVee, 
after which it falls significantly. Above 5 keVee, the total 
modulation is 0.0034±0.0024 cpd/kg, which is consistent 
with zero. The signal above 4 keVee yields a signal at 
DAMA of 0.014 ± 0.004 cpd/kg, which requires approxi- 
mately 700 kg • day of exposure for 10 events. Moreover, 
it is possible that the actual signal is at 3.5 keVee and 
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below, and the signal at apparently higher energy is due 
to the resolution of the DAMA detector [11]. 

There is a significant uncertainty in the quenching fac- 
tor as well. While q ~ 0.09 is a commonly used value, 
the measurements are uncertain, and values q = 0.10 and 
slightly higher are possible. Since the range of WIMP 
parameters allowed arises from fitting the DAMA peak, 
the uncertainty in this factor is hidden from our analy- 
ses here. Nonetheless, the presence of a larger quenching 
factor would result in a lower range of energies for the 
signal. Thus, it is clear that a robust test of the DAMA 
result involves pushing the energy threshold as low as 
possible. While the models that we consider generally 
do have signal above 50 keVr, this cannot be guaran- 
teed, especially in situations where form factors might 
suppress the higher energy events [23] . In the event the 
experiments as we describe are performed and no signal is 
seen, we would advocate lowering the threshold, even at 
the cost of exposure from reduced pressure, to whatever 
extent possible. 

With these important caveats in mind, we can proceed 
to discuss the details of what such an experiment would 
look like. 



A. Experimental Design 

Gaseous detectors can resolve the nuclear recoil tracks, 
which have lengths of several millimeters at sufficiently 
low pressures. Several gaseous directional detection ex- 
periments are already underway, including DMTPC [27| . 
NEWAGE DRIFT [H], and MIMAC fs^, which 

employ time-projection chambers to reconstruct tracks. 
However, these experiments are typically focused on spin- 
dependent WIMP-nucleus interactions and use the gas 
CF4 as a detector, with the exception of DRIFT, which 
uses CS2. For a review of the various detector technolo- 



gies, see [sTI. l32j. 

We suggest using a gas containing xenon or an- 
other heavy element. This increases sensitivity to spin- 
independent interactions because scattering rates are 
kinematically highly suppressed for lighter nuclei in the 
iDM scenario, in addition to the overall factor of 
that appears in the cross section. However, heavier ele- 
ments have shorter recoil tracks and are thus more dif- 
ficult to resolve. Furthermore, the gas should allow for 
good electron (or ion) drift and also have good scintil- 
lation properties (at least for DMTPC). Choosing a gas 
will involve some compromise between these properties. 
We note that for a splitting of 5 ~ 120 keV, A must 
be greater than 75 to see any signal for the mass range 
- 100 - 1000 GeV for an earth velocity of 225 km/s 
and an escape velocity of 500 km/s. 

According to preliminary work of the directional de- 
tection experiments mentioned above, in order to resolve 
the angles of the tracks, the gas chamber must be at a 
pressure of around 50 torr. Furthermore if the recoil en- 
ergies are too low (below ~ 50 keVr), it is difficult to 
detect the sense (head-tail discrimination) of the track, 
which reduces sensitivity significantly [H, [s^, . The 
directional resolution of the DMTPC collaboration is es- 
timated to be around 15 degrees at 100 keVr and im- 
proves by several degrees at higher energies [36| . 

The dominant irreducible background is neutron re- 
coils arising from radioactive materials near or in the 
detector. Simulations suggest background rejection is ex- 
cellent forgamma-rays, electrons, and a's [33| (see also 
Fig. 7of [23|). The DRIFT collaboration has reported on 
neutron backgrounds; however, they found a radioactive 
source {^^'^Rn) inside the detector [H. The NEWAGE 
experiment at Kamioka estimated their primary back- 
ground to come from the fast neutron flux which, when 
shielded by 50 cm of water, would contribute only a few 
events per year (39j . 



III. RECOIL SPECTRUM 



We derive the differential nuclear recoil spectrum in recoil energy En and cos 7, which is defined as cos 7 = ve -vr- 
The Earth's motion in the halo rest frame is ve and the vector vn is the nuclear recoil velocity in the Earth's frame. 
Let V be the incoming WIMP velocity in the Earth's frame. 

The single nuclcon scattering cross section is: 

da- = ^ dER dcos7 (5^^' (v - vr - (3) 

where /i„ is the WIMP-nucleon reduced mass and a„ is a reference cross section that is assumed to be the same for all 
nuclcons. The minimum velocity Vmin for a WIMP to scatter with a nuclear recoil of energy Er was given in Eq. [T] 
The differential recoil rate is for WIMP-nucleus scattering is 



= Nt— / d^v V f{v + VE)- 



dERd cos J J dERd cos j 

where /(^/), the WIMP distribution in the galaxy frame, is boosted to the Earth frame. Define the constant k: 

PX (TnTUN ifpZ + (A - Z)f, ^2 
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FIG. 2: Exposure to get a 3a measurement of (0037) 95% of the time the experiment is conducted on Earth. The energy range 
of the experiment is Eh £ [50, 80] keVr. (IRbg is the background rate; the DAMA unmodulated background rate is indicated 
by the sohd vertical line at 0.1. The bands shown give the exposures necessary as the rates modulate throughout a year. Since 
the annual modulation is asymmetric in summer and winter for low mass dark matter, the average exposure for = 70 GeV 
is indicated by the dashed line. In (a) we show three mass benchmarks from Table |I] and in (b) we show the effect of decreasing 
the angular resolution of the detector to 60 degrees and of lowering the energy threshold to 40 keVr. (Darker regions indicate 
where the bands overlap.) 



Changing variables to v' = v + ve gives: 

dR 



(i-Efld cos 7 



= {En) / d\ f{v) S^^^ {v ■ VR - VE ■ VR - f ™„ {Er, 5)) (5) 



and F^{Er) is the Helm form factor given in This formula is discussed in detail (in the context of Radon 

transforms) in [4l| . Thus we can see that at fixed Er, the signal peaks where the delta function is nonzero over 
the largest portion of the phase space, or cos 7 — ve ■ vr^ — —1. The peak in Er and fixed 7 is determined by the 
competition between the form factor (which pushes the signal to lower energies) and the inelasticity (whereby the 
minimum velocity produces a minimum value of Er). 

For a truncated Maxwell-Boltzmann distribution, the spectrum is: 

dR kF^{Er) ^( ( {vE-VR+Vrain{ER,6)f \ ( «Lc\\p,. ^ ^ cyUR\ 

- ~- ' exp 2 - exp 2~ ©(I'esc - VE - VR- Vmin{ER, 5)) (6) 



dERdcos^ n{vo) 

where n{vo) is normalization off for finite escape velocity. In the limit of (5 — > and Vesc — * oo this reduces to the 
elastic scattering result from |25|. The values we use for the astrophysical parameters are: vq — 220 km/s, ve = 225 
km/s, Vesc = 500 — 600 km/s [421, and p-^ = 0.3 GeV/cm'^. The normalized rate spectrum of several benchmark 
models is shown in Fig. [TJ 



IV. SENSITIVITY 

A robust detection of a directional modulation is pos- 
sible with surprisingly few events, and does not require 
use of the rate formulas in the previous section. In fact, a 
full likelihood analysis based on the correct model is only 
a factor of ^ 2 better than a simple technique, and for 
a convincing detection, simpler is better. In this section 
we assume the detection gas has A = 127 (for iodine; Xe 
with j4 = 131 would be similar) and focus on the energy 
range Er e [50, 80] keVr. 



A. Detectability 

For a model-independent statistic we follow [H, l43j 
and use the dipole of the recoil direction, (cos 7) . This is 
motivated by the fact that the rate should depend only 
on cos 7 and Er, so the directional part can be expanded 
in spherical harmonics. 

Our detection criterion is a measurement of (0037) 
that is 3(7 away from the distribution of (cos 7) for the 
same number of randomly distributed events. For a fixed 
exposure, we generate many random sets of model data 
(constrained by the DAMA benchmarks in Table [l|, and 
then demand that 95% of the time the result is 3a from 
the null hypothesis. The background is modeled as uni- 
form in recoil energy and angle. We assume the detector 
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has an angular resolution of 15 degrees. 

In Fig. Ufa) we show the exposures necessary for such 
conditions, as a function of the background rate, for a 
few benchmark models. At zero background, roughly 11 
events are needed for all benchmark models, on average. 
Fig. [2tb) shows the effect of decreasing the angular reso- 
lution to 60 degrees and lowering the energy threshold of 
the experiment to Ejf = 40 keVr. Because of the sharp 
angular profile of the recoil spectrum, a poor angular 
resolution does not significantly reduce the possibility of 
a detection. However, achieving an energy threshold of 
30-40 keVr dramatically lowers the necessary exposures 
because the peak of the recoil spectrum occurs at 30-40 
keVr and falls off exponentially. 



B. Parameter Estimation 



In Figs. [3][9] we show confidence levels of (68, 90, 95, 
99, and 99.9%) on the WIMP parameters for an exposure 
of 1000 kg • day. To obtain the probability, or likelihood, 
at a point in the — S plane, we either: 1) find the 
likelihood as a function of cr„ and maximize with respect 
to cr„ or 2) assume cr„ is exactly known from some other 
experiment. We can do the same also for points in — 
(T„ plane and (t„ — 5 plane. The full log likelihood function 
lives in the full 3 dimensional parameter space. Here we 
show possible slices through that space. 

For each possible slice, we have shown several vari- 
ations on the real WIMP parameters or experimental 
parameters. In the default scenario, we consider the 
m^=150 GeV benchmark with £^(^=50 keVr, a back- 
ground rate of cIRbg = cpd/kg/keVr, and v^sc = 
500 km/s. We consider the following independent varia- 
tions: 



We also perform a likelihood analysis as a measure of 
sensitivity of the experiment to the parameters of the 
model, assuming perfect energy and angular resolution. 
From our analysis in the previous section, we expect this 
does not affect the results significantly. (See also [s^ . 
which shows the sensitivity dependence on angular reso- 
lution.) The parameters we consider here are rn^,d, and 
cr„ , which we denote together simply by p. Define 



dR 



dE ad cos ^ 



{x;p)+dRBG/'2, 



(7) 



which is the rate (cpd/kg/keVr per COS7) at a given recoil 
energy and angle (denoted together by x) for parameters 
p. We assume the background rate, dRsG, in units of 
cpd/ keVr/kg, is known. 

The likelihood is the probability of parameters p given 
the events {xi}. Given events {x^}, bin the events such 
that in each bin there is only or 1 event and label the 
bins with one count by {^'i^a} and the empty bins by 
{Xfj}. The expected number of comits in a bin is 

E{X]p) ^ £li{x;p)/:^x (8) 

where £ is the exposure. Then the (log) likelihood is 



\nCtot{p) = y^ln 



a (3 

(9) 

which is the log of the Poisson probability of obtaining 
or 1 event in each bin. To find the expected average 
InCtot for a given exposure £ and true parameters poj we 
compute 



liiCtotip) =^ J dx(^fi{x;po) In fi{x;p) ~ fi{x;p) 



(10) 



which is the continuum, noiseless limit of Eq. [5] Since 
we can only compare differences in log likelihood, in this 
equation we have subtracted an arbitrary constant in p 
which takes care of the units in ln/i(a;,p). 



• Lower energy threshold (Eth — > 40 keVr) 

• Higher background {dRsG ^ 10^^ cpd/kg/keVr) 

• Higher escape velocity (uesc = 600 km/s) 

• Lower WIMP mass (m^ 70 GeV benchmark) 

• Higher WIMP mass (m^ — * 700 GeV benchmark) 

In each case, as and Vesc vary, cr„ and 5 are ad- 
justed to agree with benchmark fits to DAMA, using 
the parameters in Table [D At masses above 250 GeV, 
there is increasing tension between the DAMA result 
and other experiments, notably CDMS. This tension is 
highly dependent on the high velocity tail of the WIMP 
velocity distribution, and can be alleviated by consid- 
ering non-Maxwellian velocity distributions, for instance 
from the Via Lactea simulation 0, |45l|. Thus, we con- 
sider these points, but it should be emphasized that the 
non-Maxwellian halos generally tend to lead to a larger 
signal at DAMA (relative to the other experiments), and 
thus on a xenon target (because of the similar kinemat- 
ics), and thus we expect that our use of a Maxwellian 
distribution is conservative for these points. 

At these masses much larger than the nucleus mass, the 
threshold velocity Vmin is independent of mass and the 
spectrum depends on only through the local WIMP 
density p^/m^. In these regions and tT„ are com- 
pletely degenerate since only the combination p^anjm^ 
ever appears, as a prefactor determining the overall rate. 
This can be clearly seen in Fig.[3l which shows confidence 
intervals in the — an plane. Note that because the 
contours never close, we have have imposed the (rather 
conservative) constraint that < 100 TeV based on 
the unitarity bound |4Q] for a thermal relic. 

The effects of the — (t„ degeneracy can also be seen 
in the — 6 plane, shown in Fig. [5] Here high masses 
are all equally likely (given a fixed S) because (T„ can be 
adjusted accordingly. 

In the S—Un plane. Fig. [3 there is a sharp discontinuity 
since low masses are favored at smaller tT„ and very high 
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masses are favored at high it„. This is because at low 
scattering cross section, in order to boost the rates such 
that it matches the observed number of events, one can 
lower S or adjust the mass to optimize the number of 
rates. (The scattering rate is maximized when the mass 
of the WIMP ~ the mass of the nuclei.) However, at 
high scattering cross section, one can increase S but only 
increase the mass to very high masses to reduce the rates. 
Though lowering the mass drastically also decreases the 
rate, the angular shape at very low masses is very distinct 
(see Fig. [T]) and thus unfavored. The cutoff in Fig. [7] at 
high an is a result of the unitarity bound on the mass. 

These effects can make it difficult to constrain the 
WIMP mass at low exposures; however, it is easier to con- 
strain the ratio m^/cTn, which we have shown in Fig. [9l 

Finally, we note that in these figures we have assumed 
the earth velocity is unmodulated. For the benchmark 
where — 70 GcV, our worst-case scenario, the effects 
of the annual modulation in velocity improve the confi- 
dence levels significantly if the experiment is done during 
the summer. 

The disadvantage of the likelihood analysis is its model 
dependence. We used the truncated Maxwell-Boltzmann 
profile, whereas in reality it is likely there is more struc- 
ture in the dark matter profile. However we expect the re- 
sults to roughly be the same for many more complicated 
velocity distributions. Because of the velocity threshold 
due to 6, the inelastic scenario is almost independent of 
most streams because most streams are below the thresh- 
old velocity. Furthermore, anisotropics in the halo profile 
do not significantly affect the results here. To see the ef- 
fect of using less simplistic halo models on the elastic 
scattering spectrum and sensitivity, see [i^l and [ist . 

V. CONCLUSIONS 

Motivated by the finding [20|] that inelastic dark matter 
(iDM) is compatible with both the DAMA annual modu- 
lation signal at 22 — 66 keVr and limits from other exper- 
iments at lower energies, we have investigated prospects 
for directional detection in the context of the iDM model. 
We are encouraged by the fact that ZEPLIN-III has 
also detected a number of events in the 40 — 80 kcVr 
range 49]. This has not been claimed as evidence of 
WIMP scattering, but makes it impossible to rule out 
iDM with such data. In the near future, LUX (5^ and 
XENONIOO [50] will have greatly improved sensitivity 
and lower backgrounds, and will provide a sharp test of 
the iDM/DAMA scenario. If these experiments also de- 
tect an excess of events above background in the appro- 



priate energy range, a major effort in directional detec- 
tion will be justified. 

Directional detection with a gaseous detector contain- 
ing a heavy gas (e.g. Xe) may not require the huge 
exposure times implied by the elastic scattering limits. 
For a threshold energy of Eth = 50 keVr, we find that 
exposures of order 1000 kg • day in a directional ex- 
periment can convincingly refute or support the claims of 
DAMA in the context of the inelastic dark matter model. 
At zero background, roughly 10 events are needed for a 
clear detection of WIMP scattering. Even with larger 
backgrounds, the required exposure is a few hundred 
kg • day, over most of the iDM parameter space that 
can explain both DAMA and other direct detection ex- 
periments. With roughly 1000 kg • day, it is possible to 
obtain a measurement of (5 > at high significance and 
also the parameter m^/cr„ to within an order of magni- 
tude. 

Furthermore, if it is possible to roughly determine one 
of the WIMP parameters, for example 6 ~ 120 keV, 
via another experiment, the mass and nucleon scatter- 
ing cross section are highly constrained with an exposure 
of a few hundred kg • day because of the distinctive shape 
of the energy-angle recoil spectrum. 

Significantly lower exposures are needed if the thresh- 
old energy is decreased. As discussed in Section |TT1 be- 
cause of the uncertainties in the nuclear recoil energies 
of the DAMA signal, it is crucial to reduce the thresh- 
old energy as much as possible. For low masses, the re- 
coil spectrum is sharply distributed in energy and angle. 
However, typical recoil energies are smaller. Thus with 
an energy threshold of Eth = 50 keVr most of the events 
for ~ 70 GeV are not seen. With an energy threshold 
of 100 keVr and = 70 GeV, none of the WIMP re- 
coils can be seen. Though the required volume increases 
and angular resolution decreases when Eth is lowered, 
we found that a poor angular resolution (~ 60°) does 
not significantly affect the results, assuming that 3D re- 
construction of the track and determining the sense is 
still possible. 

We acknowledge helpful discussions about directional 
detection with Peter Fisher, Jocelyn Monroe, and 
Gabriella ScioUa. Rick Gaitskell, Dan McKinsey, and 
Peter Sorensen provided helpful advice and much-needed 
skepticism. DPF is partially supported by NASA LTSA 
grant NAG5-12972. TL is supported by an NSF Grad- 
uate Fellowship. This work was partially supported by 
the Director, Office of Science, of the U.S. Department of 
Energy under Contract No. DE-AC02-05CH11231. NW 
is supported by NSF CAREER grant PHY-0449818 and 
DOE OJI grant #DE-FG02-06E R41417. 



[1] R. Bernabei et al. (DAMA), Phys. Lett. B480, 23 (2000). 

[2] R. J. Gaitskell, Ann. Rev. Nucl' Part. Sci. 54, 315 (2004). 

[3] J. Angle et al. (XENONIO), Phys. Rev. Lett. 100, 021303 (2008), 0706.0039. 

[4] J. Angle et al., Phys. Rev. Lett. 101, 091301 (2008), 0805.2939. 



8 






-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

l0g,o (On/Oo) 



-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 



-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 






-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 
log,o {oJo„) 



-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 
log,o (oJa„) 



-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 
log,o (On/Oo) 



FIG. 3; Confidence levels for determining nix s-nd ffn, where S is unknown, with an exposure of 1000 kg • day. ao = 10 '^'^crn^ 
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FIG. 4: Confidence levels for determining and cr„, where 5 is known with an exposure of 1000 kg • day. ctq = 10 *°cm^. 






FIG. 6: Confidence levels for determining and 5, where cr„ is known, with an exposure of 1000 kg • day. 
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FIG. 7: Confidence levels for determining 5 and a, where nix is unknown, with an exposure of 1000 kg • day. ao = 10 ^'^cm^ 
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FIG. 9: Confidence levels for determining S and m^/an, where rn^ is unknown, with an exposure of 1000 kg • day, taking 
(Jo = 10^*" cm'^. Over most of tlie parameter space, some value of rn^ (and therefore an) can be found to produce enough 
events for the given 6. However, in the case of large S and large mx/crn, no solution is possible in some CEises. 
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